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Abstract 

The existing equivalence between (generalized) BMS Ward iden¬ 
tities with leading and subleading soft graviton theorems is extended 
to the case where the scattering particles are massive scalars. By ex¬ 
tending the action of generalized BMS group off null infinity at late 
times, we show that there is a natural action of such group not only 
on the radiative data at null infinity but also on the scattering data of 
the massive scalar field. This leads to a formulation of Ward identities 
associated to the generalized BMS group when the scattering states 
are massive scalars or massless gravitons and we show that these Ward 
identities are equivalent to the leading and subleading soft graviton 
theorems. 


1 Introduction 

Since the seminal work of Strominger DEa. there has been a resurgence of 
interest in analyzing the role played by BMS group |3] and its various ex¬ 
tensions |1] as symmetry group of (perturbative) Quantum Gravity S-matrix 
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[HE]. However work done so far is restricted to pure gravity or pure grav¬ 
ity coupled to massless matter. In this paper, we extend the relationships 
between a group of asymptotic symmetries Q defined in [TJ [8] (referred to as 
generalized BMS group) with leading and sub-leading soft graviton theorems 

PHD]. 

Our work is based on the ideas presented in HD, where the Weinberg soft 
photon theorem for scalar QED with massive charged particles was shown 
to be equivalent to Ward identities associated to large U{1) gauge transfor¬ 
mations which acted on the asymptotic phase space of massive scalar field 
defined at time-like infinity. In this paper we show that there exists an ac¬ 
tion of Q on the asymptotic phase space of massive scalars such that charges 
associated to Q are sum of charges associated to time-like inhnity and null 
inhnityj^ We then show that the statement of Ward identities corresponding 
to such charges is equivalent to leading as well as sub-leading soft graviton 
theorems. 

A key player in obtaining this equivalence are certain ‘boundary to bulk’ 
Green’s functions that map generators of Q (which are vector fields at null 
infinity) to the asymptotic (at time-like infinity) bulk vector fields. We show 
that these Green’s functions are intricately tied to the soft factors which 
arise in the soft theorems and this fact plays the central role in obtaining 
the equivalence between soft theorems and Ward identities. Further details 
of these Green’s functions appear in the companion note [H]. 

The outline of this paper is as follows. After giving a brief conceptual sketch 
of the key ideas which underlie the proof (establishing equivalence between 
Ward identities and soft theorems), in section]^ we establish the asymptotic 
phase space T of the scalar field-gravity system. This phase space is a direct 
product of radiative phase space of gravity defined at null infinity and the 
asymptotic phase space of the massive scalar field defined at time-like infin¬ 
ity. In section]^ we show that there is a group of large diffeomorphisms which 
are non-trivial at time-like infinity and are (i) obtained from the generalized 
BMS vector fields localized at null-infinity with help of certain Green’s func¬ 
tion with well-defined boundary conditions and (ii) preserve the de Bonder 
gauge-fixing condition of perturbative gravity. This gives us a unified picture 
of the asymptotic symmetry group of gravity associated to null and time-like 
inhnity. In section it is shown that the action of this symmetry group on 

^For the supertanslation subgroup ST C G, the action on the massive scalar seems to 
agree with the one proposed by Longhi and Materassi m (see also m)- 
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the total phase space T is symplectic and we derive the charges associated to 
the symmetry generators. In sections and we establish the equivalence 
between soft theorems and Ward identities associated to such charges. We 
end with some remarks and conclusions. 

1.1 Basic sketch of the proof underlying the equiva¬ 
lence 

We now illustrate the sequence of steps which leads to the equivalence in 
the context of the supertranslation subgroup ST C ^ and Weinberg’s soft 
graviton theorem. It is important to recall that the underlying theory is 
perturbative quantum gravity coupled to a massive scalar held where we use 
de Bonder gauge to describe the linearized metric. 

Let us for concreteness focus on future asymptotics. We use two different 
systems of coordinates: (n, r, x) adapted to future null inhnity and (r, p, x) 
adapted to future time-like inhnity Null inhnity is reached by taking 
r —)■ cx) in the hrst coordinates and time-like inhnity is reached by taking 
r —)■ cx) in the second coordinate system. 

Given a supertranslation vector held fdu at X'*', we can consider an as¬ 
sociated spacetime ‘bulk’ vector held that is a residual ‘gauge’ symmetry 
in de Bonder gauge, i.e. satishes = 0 and its boundary value at 

X’*' is fdu- At time-like inhnity such vector held has the asymptotic form 

~ leading component can be determined via a 

Green’s function G{p,x;y) that maps f{y) to by: 

tf{p,x)= j £yG{p,x;y)f{y). ( 1 ) 

J 

If we now consider a scattering amplitude involving n massive particles 
then the corresponding Ward identity for supertranslation [2], 

(out|[Q®°^*,5]|in) = — (out I 5] I in) (2) 

takes the form (for details see the main text of the paper): 

E f 

lim — / (f^w f(w,w)D‘^ {out\a+ (Eg, w,w)S lin) = 

Es^0+ ZTl J 

n 

i=l 


j (^yG{\pi/mi\,pi;y) f{y){out\S\m). ( 3 ) 
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It is important to notice here that in the r —?• oo limit the momentum of a 
free massive particle p determines the point on T-L^ the particle reaches via 
p = \p/m\, X = p. 

This formula looks structurally similar to Weinberg’s soft graviton theo¬ 
rem which is given by 

^ + \2 

lim Es{ont\a+{Es,w,w)S\m) = y^ . „ f' —(out|5|in) (4) 

Es^o+ ^[q/Es)-pi 


We show below that after expressing boundary coordinate y into sphere co¬ 
ordinates {w,w), one has the following remarkable relation between Green’s 
function associated to supertranslation generators and soft factor 

Pr? 


1 

-miG{\pilmi\,Pi]w,w) = —Dl - 

27r [q/Es) ■ Pi 


(5) 


This relation is the key in establishing the equivalence between ST Ward 
identities and Weinberg’s theorem. Similar sequence of logic interspersed 
with Green’s function associated to sphere vector helds on the boundary 
and certain derivatives of sub-leading soft factor leads to the equivalence 
between Ward identities associated to Diff(S'^) C Q and Gachazo-Strominger 
(GS) subleading soft theorem. 

We conclude by commenting on the relationship with the case of massless 
particles. The equivalence in such case can be cast in the above language 
by expressing the supertranslation action on massless particles in terms of a 
‘Green’s function’ that is just an identity kernel. The analogue of relation 
(|^ is then: 

- Ei5^‘^\zi,w) =—Dl, ' — (massless particles), (6) 

27r [q/Es) ■ Pi 


where (zj, Zi) are sphere coordinates for p*. The p —)■ oo boundary behavior 
obeyed by the Green’s function G(p, £; y) ensures that Eq. (|^ reduces to 
Eq. ([^ in the m —)■ 0, p = constant limit. 


2 Asymptotic phase space 

We consider perturbative gravity coupled to a massive scalar held in de Bon¬ 
der gauge. We assume a total phase space of the form T = Tgrav x T^ with 
Tgrav and T^ the free-held asymptotic phase spaces of gravity and massive 
held respectively. 
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2.1 Gravity phase space Fg^av 

The gravity phase space is described by the ‘radiative data’ Cab{u, x) at null 
inhnity with symplectic structur^ 

f^grav(5 ,= [sC^^S'Cab -S^d'Y (7) 

In quantum theory this data is related to the asymptotic Fock functions as 
follows. Consider the Fourier transform of Cab, 

poo 

Cab{E,x) := / CAB{u,x)C^^du, ( 8 ) 


and go to {z, z) coordinates on the sphere. Then for > 0 the positive and 
negative helicity graviton anihilation functions of momentum p = Ex are: 


27r'Z 27rz 

a+{E,x) =—C^^{E,x), a-{E,x) =—C 2 z{E,x), (9) 

where = 2/(1 + zzY- The symplectic structure ([^ implies the standard 
linearized gravity Poisson brackets: 


{ah{^, = -i2Ep5hh'{'^'xfd^^\p- P)- ( 10 ) 


2.2 Massive scalar field phase space F^^ 

The massive scalar field phase space is described in terms of data on a unit 
hyperboloid "H describing time-like infinity. The coordinates adapted to such 
description are 

r := yjE — C, p := ^ , (11) 

in terms of which the Minkowski metric reads: 

-I- (12) 

with 

do^ 

da"^ = -r -I- p^'jABdx^dx^ =: haBdx°‘dx^ (13) 

I + 

"^CABiu.x) is assumed to satisfy Cab{u,x) = C^g{x) + 0 (| m |“'^) as u ^ ±oo. 
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the unit hyperboloid metric (with scalar curvature -6). At large r the free 
massive held behaves as: 


p, i) = (i>(P, + i>'(P. *)«"”■) + (14) 

with b{p, x) representing free data. To ensure well-dehnedness of upcoming 
expressions involving integrals in we will assume the free data satishes 
‘hnite energy’ p —)■ cx) fall-offs: 

6(p,x) = (15) 


The symplectic structure is 




im 


2(27r) 


d^V {Sbd'b* - 6 ^ 6'). 


(16) 


'n 


In quantum theory b{p, x) become the anihilator operator of a scalar particle 
with momentum p = px. The symplectic structure (16) implies the standard 
Poisson brackets: 


{6(^,6*(/)} = -i{2Trf{2Ep)6^{p-:^). (17) 

3 Extension of Q to time-like infinity 

The residual gauge transformations in de Bonder gauge are generated by 
vector helds satisfying the wave equation (with respect to the hxed reference 
Minkowski metric): 

□r = 0. (18) 

We are interested in ‘large’ gauge transformation that are non-trivial at in- 
hnity. At null inhnity we would like to have generalized BMS vector helds. 
These vector helds are dehned by the condition of being asymptotically 
divergence-free at null inhnity [7]: 

VaC = O(r-i). (19) 

These vector helds are parametrized by sphere functions / (supertransla¬ 

tions) and sphere vector helds (generalized rotations) according to, 

dr, u, x) = fdu V^dA + uadu - radr + ..., (20) 
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where 2a is the 2-d divergence of and the dots indicate snbleading terms 
in the 1/r expansion. In the usual treatment these terms are determined by 
the Bondi gauge condition: 


9rr = 9rA = 0, det gAB = r* det 


( 21 ) 


where qab is the unit round sphere metric (in spherical coordinates, det qab = 
sin^6*). However in the present case one should instead use condition (18) 
to £x such subleading term (see [12] for such determination in the case of 
supertranslations). For the purposes of the present paper we will not need 
the specihc form of such sub leading terms. 

We start with the following ansatz for the asymptotic expansion of the 
vector helds off time-like inhnity 


C(t,p,x) = r(P.4)+0(T') 

5“(t,P, 4) = r(p,i) + 0(T->) 


( 22 ) 

(23) 


where a, (3, ... denote indices on the hyperboloid. As we will see shortly, this 


ansatz is consistent with Eqs. ( [l^ , ( [19| ), ( |^ . Note also that a vector held 
satisfying 


has a well-dehned action on the massive held free data 
b{p,x) (obtained by evaluating the derivative ^°'daP> in the expansion ([l4|) 
given by: 

S^b = -imCb + Cdab. (24) 


The idea is to use (18), (19), (20) to determine and in terms of / 
and V^. In appendix |A| it is shown that conditions (18), (19), (20) imply: 


Ar = sr, 


limp C{p,x) = f{x), 

p^OO 


(25) 


Ar = 2r. C„r = 0, limp(9,i) = l/-i(i), (26) 

P^OO 

where A and are the Laplacian and covariant derivative on 'H. If G(p, g) 


and GJKp, f; g) are Green’s functions for equations (25) and ( |26| ) respectively, 
the solutions can be written as: 


Cip,x)= d^qG{p,X]q)f{q) =-. fn{p,x), (27) 

Js^ 


r(p,£)=/ c/^gG^(p,a;;g)H^(g)=:H“(p,x). (28) 

Js^ 
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The explicit form of G{p,x-,q) and f;g) is given below in Eqns. (73) 
and (106) respectively. Fnrther details of these Green’s fnnctions are given 

in [H]. 

The above extension of generalized BMS vector helds to time-like inhnity 
provides an action of generalized BMS on the phase space of massive particles 
and hence on the total phase space T = Tgrav x T^- We now describe the 
associated charges. 


4 Action of Q on P 

In this section, we derive the complete expression for charges associated to 
all the generators of Q. As T is the Cartesian prodnct of radiative phase 
space of gravity dehned at I and asymptotic phase space of scalar field 
dehned at "H, charges associated to any (generalized) BMS vector held are a 
snm of gravitational charges dehned at X ElE], and scalar charges dehned 
at Ti. Exactly as in the case of pnre gravity [3E], charges associated to 
snpertranslation and sphere vector helds split into contribntions which can 
be classihed as hard and soft charges. The soft charge is linear in an inhnite 
wavelength mode of Cab- As the gravitational charges associated to the 
radiative phase space Tgrav were derived in earlier papers, we do not reprodnce 
the derivation here and only derive the charges associated to scalar held which 
are fnnctions on T^. 


4.1 Supertranslation charges 

From the previons consideration, we have an action of snpertranslations on 
the total phase space T = Tgrav x r^: 


SfCAB = fd^CAB-iiDADsff’' 
Sjb = -imfyb, 


(29) 

(30) 


with /-^ dehned in (27). The generator Qj satisfying 5Qf = (2(5,5/) can be 
written as 

Q/ = Q5-d + Q-ft_ ^3^) 

The hard part is a snm of gravitational and matter contribntions 


^hard _ /^harcl\ i /^hard\ 

Wf — “r f Jmatter 


•^hard 


hard\ 


( 32 ) 





The gravitational contribution was derived in na and is given by 


(Q/"’'")grav = i,n,UfduC,C) 


2 

1 


= - du^fduC^^d^CAB. 


The matter contribution is given by: 


(Q5”")matter = ^^m{6fb,b) 


1 


m 


2(27r) 


d^VfuVb. 


(33) 

(34) 

(35) 

(36) 


'v. 


One can verify that (36) coincides with the expected expression in terms of 
the energy-momentum tensor: 


(Qj“'^)matter = “ hm 


dSa^T^i). 


'Hr 


(37) 


where T-Lr is the r =constant hypersurface and = /-h^t + 0{t~^). For 
the translation subgroup of supertranslations, Eq. (37) gives the total linear 
momentum of the held (in the asymptotic future). 

We now describe the soft part in (31). This was derived in [I^, and can 
be written as: 


(38) 


Qr = Ograv(-2/lV,C') = -- / d^VfD^D^^lCAB]. 

^ Js^ 

Where the square bracket denotes difference of boundary values at m = ±oo: 

[CabKx) ■= Cab{oo,x) - Cab{-oo,x). (39) 

In order to establish the equivalence between ST Ward identities and Wein¬ 
berg’s soft theorem we need to impose the additional condition]^ 




(40) 


This condition is satished by the so-called Christodolou-Klainermann space- 
times considered by Strominger in [1]. 


’Conversely, Weinberg’s soft theorem can be seen to imply condition (40). 
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We conclude the section by writing the soft charge in terms of the mode 
functions. First, we express [Cab] as a zero energy limit of the Fourier 
transform Cab{E,x): 


[CABix)] = -i lim ECab{E,x). 
E^0+ 


(41) 


Using (§ and ([40|) the soft charge ([M| can be written as: 


= — \im E / d^zfDia+{E,z,z) 

ATI E-5-0+ / 


(42) 


(or alternative expression in terms of D^a-). A technical but important 


point is that when performing covariant derivatives as in (42) one needs to 


remember the two-dimensional tensor structure (including density weight) of 
the quantity being derived. For the case of a+ this tensor structure can be 
read off from Eq. 0. Doing so one hnds that the differential operator in 


(42) acts as: 


Dja+{E,z,z) = 5j(7“5j(7.ia+(B,z,z))). 


(43) 


4.2 Sphere vector field charges 

In this section we derive the charges associated to the generators of Diff(S'^) C 
Q. The corresponding generators of Q are vector helds at null inhnity which 
are in turn parametrized by vector helds V^Oa on the conformal S"^. 

^'^{u,x)da= V^dA + uadu, (44) 

a = ^(DaV^). 

For these vector helds their action on F is given by, 

SvCab = S^^'^'^Cab + S^^Cab (45) 

dyb = Cv^b, (46) 

where 

dy^'^^CAB '■= EvCab — O'Cab + (ouOuCab (47) 

d^f^CAB := -2u{DADBaf^ (48) 


and Vbl dehned in (28|) 
speaking the action 


There is a subtlety in the present case in that strictly 
(|45) is between radiative phase spaces associated with 
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different 2-dimensional metrics Qab- As shown in [S], one can nevertheless 
compute the associated charge by embedding Fgrav into a larger space that 
allows for variation of qab- Doing so one hnds the charge is a sum of ‘hard’ 
and ‘soft’ pieces. The ‘hard’ piece turns out to coincide with the naive 
expression, 

(Q^“")grav = (49) 

= \j du^/^ duC'^^{CvCAB - oiCab + auduCAB)- (50) 

On the other hand, the soft charge receives contributions from the variation 
of the 2-metric SvQab resulting inj^ 

Qft = 1 y du^ (51) 


which differs from the naive expression f2gi.av(^y^*C', C). 

For the matter contribution, we note that the fact that is divergence- 
free as shown in Eq.(26) implies the action (46) is symplectic on Fm]^ The 


associated charge is then given by 




-QrniSyb, b) 


im 


2(27r) 


'n 


d?Vb*Cv^b. 


(52) 

(53) 


One can verify (53) coincides with the expected expression from the energy- 


momentum tensor perspective. In the notation of Eq. (37): 


= - lim / dSaT^ev. 

JUr 


(54) 


where ^ydb = + 0{t~^). For the vector helds associated to rotations 

and boosts, Eq. (|3^ gives the total angular momentum of the held (in the 
asymptotic future). 


'‘This charge is only defined on the subspace of Fgrav satisfying the stronger fall-offs 

was first given in [6] for the case where 


±oo. The charge 


51 


C{u,x) = 0(|u|“‘“'') at u 

is local conformally Killing. In such case SyqAB = 0 and the charge can be derived 
within Fgrav 

®Here we are crucially using the fact that the field is a scalar. For fields with nonzero 
spin the symplectic form will depend on the full hyperboloid metric (not just its 
volume element Vh) and further subtleties will arise. 
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The total sphere vector held charge is then given by 


Qv — Qv 

')hard'\ 


+ Q 


V 


( 55 ) 


with the hard piece = (Qy‘^^'^)grav + given by Eqns. (50), 

(53) and the soft piece given by Eq. (ig. We conclude the section by 
writing the soft charge in terms of the mode functions. Using ([^ and the 
prescription given in [6] that projects out the Weinberg pole, f duC(u,x) = 
limE_,.o(l + EdE)C{E,x), the soft charge can be written as: 


yOSOft _ 

Vv — 


— jim^(l + £aE) J i(‘z{V^Dla+(E,z,z) + V‘Dla_{E,z,z)). 

(56) 

Taking into account the tensorial structure of the mode functions, the explicit 

(67) 


action of the differential operators in (56) is found to be given by: 

L>?a+ = a|a+, Dla_ = d^a.. 


5 Supertranslation Ward Identity = Wein¬ 
berg soft theorem 

Having derived the action of generalized BMS group on the gravity-massive 
scalar held asymptotic phase space, we now turn to the quantum theory 
and analyze the constraints such a symmetry imposes on the perturbative 
S-matrix of the theory. That is, we ask that if indeed the generalized BMS 
group was a symmetry group of the perturbative S-matrix, what would be 
its implications. 

As in the case of gravity coupled to massless particles [5], in this section 
we show that the inhnity of Ward identities associated to supertranslation 
subgroup are equivalent to the Weinberg’s soft graviton theorem. We hrst 
review the equivalence of supertranslation Ward identities with the Weinberg 
soft theorem in the case of pure gravity and then generalize it to our case. 

5.1 Review for the case of external massless particles 

To setup notation for the later section we review the equivalence for the 
case of external massless particles. For given ‘in’ and ‘out’ states composed 
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of massless particles of momenta {p*}, the Ward identity associated to a 
supertranslation /, 

(out|Q/-5 — 5Q/-|in) = 0 (58) 

can be written as 


lim — / (fwf(w,w)D‘^{out\a+(Es,w,w)S\m) = 

Es^0+ ZTT / 


^Ei/(pi)(out|5|in). (59) 


In ( [5^ the sum is over all external particles, with Ei = ±|pi| for outgo¬ 
ing/incoming particles and pi = Pi/\'Pi\- 

On the other hand, Weinberg’s soft graviton theorem can be written asj^ 

lim Es(out|a+(ii^s, ta, ta)S'|in) = ^(out I S'I in), (60) 

^ {q/Es)-pi 

where q/Eg = (1,^) with q parametrized by {w,w) and e~^{w,w) the polar¬ 
ization vector [5]: 

= 1/V2{w,l, —i, —w). (61) 

If we parametrize an external momentum p by {E,z,z), the soft factor in 
(60) take the form: 


w) ■ pY 
{q/Es) ■ p 


= —E s( 2 ;, z] w, w), 


with 


s( 2 ;, w, w) := 


1 -f ww w — z 
1 -I- w — z 


( 62 ) 


(63) 


To go from ( |6b| ) (soft theorem) to (59) (Ward identity for supertranslation /) 
one performs the operation operation (27r)“^ J cPwf{w,w)D^ on both sides 
of (60). The LHS becomes the left term in (59). That the RHS also coincides 
follows from the identity (see appendix [B|): 

D'^s{z, z;w,w) = 27 i5^‘^\w — z). (64) 


®Here and below, the soft theorems are written with sign convention such that all 
particles are outgoing. 
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To go from (59) to (60) we look at the Ward identity for the particular 
function 

f{z,z)=s{z,z;w,tu). (65) 


The RHS becomes the right term in (|60|). That the LHS also coincides follows 

( 66 ) 


from an integration by parts and the identity (see appendix [B|): 

DI s(;s, z; w, w) = — z). 


5.2 External massive particles 

By using normal ordered prescription to dehne the quantum charge 
one has, 

[Kp), = mfn{p/m)b{p) (67) 

(see Eq. (30)). Using (67), the proposed Ward identity (58) for external 
massive scalars takes the form: 


lim ^ / (fwf{w,w)D‘^{out\a+{Es,w,w)S\m) = 

Es^0+ ZTl J 

-Yl'^if‘H{Pi/m){out\S\m), ( 68 ) 

i 

where rrii = ±m for outgoing/incoming particles. On the other hand, Wein¬ 
berg’s soft graviton theorem takes the same form as in Eq. (60): 

lim £'s(out|a+(E's, tc, t(;)5|in) = -— / ’ —^^(outliSlin). (69) 

Es^o+ ^ [q/Es)-pi 

We now repeat the steps that led to the equivalence in the massless case. 
If we perfom the operation ( 27 r)“^ J cPwf{w,w)D'^ on both sides of (69) we 
obtain: 


lim ^ 

Es^0+ 271 


with 


d‘^wf{w,w)D^{ont\a-^-{Es,w,w)S\m) = — mi/(pi/m)(out|5|in), 

(70) 


f{p/m) = J (fwG{p/m;w,w)f{w,w), 
Gyjpim-.w.w) ■.=-—D^- 


271 ^ {q/Es) ■ {p/m) 


(71) 

(72) 
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Parametrizing the 3-momentum particle as p = mpx one can verify that (72) 
takes the following simple form: 


G{p/m;q) = 




47r {{q/Es) ■ {p/m)) 


3 • 


(73) 


Now, by direct computation it can be verihed (73) satishes (see also 

(A(p,x) - 3)G(p, X- q) = 0. (74) 

Furthermore, using the fact that the p —?• oo, m = constant limit can be 
written as a m —)■ 0, p =constant limit, together with Eq. (64) one Ends: 


lim p ^G{p,x]q) = 5^‘^\x,q). 

p—^OO 


(75) 


It then follows that G is the Green’s function for Eq. (25), and so / = fu- 
Thus the identity (70) coincides with the Ward identity (68). 

To go from (68) to (69) we repeat the steps as in the massless case: 
Consider the Ward identity for the particular function 

f{w,w)=s{w,w]Zs,Zs). (76) 


where s{w,w-, Zs,Zs) is defined in Eq.(63). By the same argument as in the 
pure gravity case, the LHS then becomes the left term in (69). The multi¬ 
plicative term on the RHS is given by 

-,nMpM = -m/d=«,G(p7™;«..®)s(«,.®;(77) 

_ __ {e+(w,w)-(p/m)f 
ialE.) ■ (p/m) ’ 


where we used the definition (72) of G, integrated by parts and used Eq. 
(66). Thus one recovers the right term in (69). 


6 Sphere vector field Ward identity = CS soft 
theorem 

We now turn our attention to the relationship between Ward identities arising 
from the Diff(S'^) vector helds and the generalization of Cachazo Strominger 
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soft theorem to the case where the external particles are massive scalars. Such 
a theorem is a factorization formula which relates the scattering amplitude 
involving scattering states that include massive scalars as well as one sub¬ 
leading soft graviton in terms of certain soft factors and amplitudes involving 
only the external massive scalar particles. This formula has only been derived 
in the limit that masses of scattering particles is zero. However following [12] 
we will try to give a persuasive argument to the effect that such sub-leading 
theorem continues to hold in the present case as well. 

After presenting such argument we will review the Ward identity-soft the¬ 
orem equivalence for massless particles and hnally establish the equivalence 
for massive particles. 


6.1 Subleading soft theorem for external massive scalars 


The set up is as follows. We consider scattering amplitudes containing ex¬ 
ternal scalar particles and one graviton and are interested in the limit of such 

scattering amplitudes given by hm£;^_> o -Mn+iiPi, ■ ■ ■ ,Pn, Qs) — 


Here is the Weinberg soft factor whose structural form re¬ 

mains the same irrespective of whether pds are time-like or null. Sub leading 
soft theorems like CS theorem are statements which express the above quan¬ 
tity in terms of Ain in fhe sense that 


lim 

0 


Mn+l{pi,.. . ,PnAs) 


Eg 


S^^^Ain. (79) 


In the case that the scattering particles are massless, in a beautiful pa¬ 
per Broedel et. al. [T6| prove that by postulating the form of as 

a differential operators on the momentum space of scattering particles]^ 
one could constraint their form severely by demanding that scattering am¬ 
plitude to be Poincare invariant, gauge invariant and that these differen¬ 
tial operators associated to split as a sum over external particles as 
^( 1 ) 

= 9pJ. The last condition is what is referred to as 

the locality constraint and is a reasonable condition to require when we are 
working with tree level amplitudes. As noticed in [T2] there is yet another 
non-trivial constraint the soft factors have to satisfy which is related to the 

^Their arguments in fact are more general and apply to scattering particles with spin as 
well (in which case the postulated form of 5^^ involves differential operators in momentum 
space as well as differential operators on the space of polarization tensors). 
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fact that scattering amplitudes are distributions in the momentum space as 
there is always an overall momentum conserving ^-distribution multiplicative 
factor in their dehnition. If we refer to the scattering amplitudes without the 
momentum conserving (5-fn. factor as reduced scattering amplitude, then 
by assuming that the soft factors are the same for the reduced as well as 
un-reduced scattering amplitude^ it was shown that has to satisfy, 

■ dp6^^\P)). (80) 

The most general form for = '^iS^^\e{q),q,Pi,dp^) can be postulated 
based on the following conditions: 

(i) Poincare invariance implies that has to be linear in polarization ten¬ 
sor tensor, 

(ii) Dimensional analysis shows that must have mass dimension zero 
and 

(iii) As it is sub-leading it must be invariant under g —?• X q for scaling 
parameter A. 


Condition (i) in conjunction with locality requirement implies that 

s'” = (81) 

where is a multiplicative factor which doesnot involve differential 

operators on momentum space. 

If the scattering particles are massless, then the most general form for 
jg consistent with conditions (ii), (iii) above is [TB] 




E. 


1 (g.pi)2 


+ c! 


' (Q-Pi) 


+ d 


rj'^^P-pX q° 
q-Pi) 




(82) 

As shown in [16], gauge invariance and distributional constraints imply 
that in the massless case (that is when Pi-pi = 0), is zero, the higher 

derivative operators are absent as distributional constraint implies that we 
require each of them to annihilate the delta-function and the remaining hrst 


®This assumption was not needed for Broedel et. al. and this fact is actually a con¬ 
sequence of all the constraints we have listed above. However this fact is not relevant for 
the issues being addressed here and hence we ignore it for simplicity. 
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order differential term precisely reduces to the Cachazo-Strominger soft fac¬ 
tor. We will now show that this result remains uneffected even when we drop 
the Pi - Pi = 0 condition. 

As in massless case we postulate that the factorized sub-leading soft factor 
takes the form 




function 


(83) 


where is multiplicative and does not involve differential operators. 

The hrst time involves linear derivative operators as well as operators that 
will be 0{d‘f). However exactly as in the massless case, the distributional 
constraint will imply that these higher order operators have to annihilate 
5^{YliiPi) and hence will be zero. Thus just as in the massless case, the func¬ 
tional form of S\^\e^y^q^pi^di) is given by the same expression that appear 
in Eq.rtsTj). 


As the vectors Pi are not null and if we assume that all the scattered parti¬ 
cles have same mass m then the most general form of consistent with 

conditions (ii), (iii) stated above, and such that on contraction with and 
p1 ^ no two terms get repeated is given by 


u 


fiupa _ 




+ 


cirri'^ - 




(pi-qP 


+ c; 


6 {pi-q)rrfl 


' ^7 m2 


(84) 


where ... indicate the terms inside the parenthesis in eq.(82). 

We now impose conditions of gauge invariance and show that all the 
terms in Eq.(84) which were not present in Eq.(82) necessarily vanish. We 
immediately see that owing to traceless-ness of the term associated to c\ 
is absent. 

Gauge invariance of scattering amplitude implies that we need to impose 

= A(^g,)) = 0 (85) 

for any gauge parameter A^. This condition gives us 


E, (2cl + + ci)qPq‘^^ + {cl +ci)qPX‘^ + {cl + 4)g-A^ 


q-Pi 


I (Pi-A) p a I i + {pi-q)pPy^ 

-r Cq ^2 q Pi -r Cy ^2 


( 86 ) 


nr-^M — 0 
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We thus have following set of equations which includes the set that was 
obtained in [16] 

2c\ + C 2 + Cg = 0 

4+ 4 = -(4 + 4) = c 


4 = 0 

C*7 = 0 


(87) 


The second equation allows that some of the terms in eq.(86) vanish due to 
conservation of total angular momentum. 

Gauge-invariance also shows that Sr^Ltion vanishes. As we require 


4 

is 


( 1 ) 


function 


to satisfy conditions (i),(ii),(iii) listed above, its most general form 


c(l) - V 4 

"^function ~ ^8 


U u 

Pi Pi 


9 


Clearly gauge invariance would imply that 


( 88 ) 


(89) 


which for arbitrary A, imply that Cg = 0. 

Thus, the extra possible terms with coefficients c^, Cq, cy, Cg all vanish and 
the soft factor takes the same form as in the massless case. That is, due 
to Poincare invariance, locality, gauge invariance as well as distributional 
constraint, assuming that there is a factorization in the sub-leading soft limit, 
the corresponding sub-leading soft factor takes the same functional form 
irrespective of the masses of the hard particles. 


6.2 Review for the case of external gravitons 

We first sketch the equivalence between sphere vector field Ward identities 
and CS soft theorem for the case of external gravitons. We refer to [7| for 
details. The Ward identity that follows from the condition 

(out|Qy5 — 5Qy|in) = 0 (90) 
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for and ‘in’ and ‘out’ states composed of gravitons of momenta {pi} can be 
written as: 


lim {l + EsdE,) 
dvr Es^o+ 


d^w{V'^d:^{out\a+{Es,w,w)S\m) + f/"’5;^(out|a_(ii^s, it','W^)‘5|in)) = 

Jy(out|5|in). (91) 


Here Jy is a differential operator associated to that acts on the momen¬ 
tum variables of the z-th graviton. Its form is obtained from the action (47) 
on the mode functions ([^ (see Eq. (57) of [7] for the explicit expression). 

On the other hand, Cachazo-Strominger subleading soft theorem can be 
written as: 

lim (l-Fi7sClEj(out|a)|)^*(i7s,r(;,-u})iS|in) = £+(«;,M})gi,Jf'^(out|5|in), 

! _i.n+ f ^ ^ ^ 


Es^0+ 


Q-Pi 


(92) 

with and q as in Eq. (60) and with the total angular momentum 
operator on the z-th graviton |j In [7] we noticed that the soft factor in (92) 
can be written as: 


e^(w,w)■p^ 


q-Pi 




KT -P 


(93) 


where J* + is the differential operator associated to the sphere vector field 

{w,w) 


K. 


w 


{w,w) 


-ds. 


w] 


(94) 


Similarly the factor associated to a negative helicity subleading soft graviton 
can be written as the differential operator associated to the vector held 
that is the complex conjugated of (94). 


®We follow the conventions used in m where — Pfj.d^+ helicity terms. 

There are sign errors in Eqns. (63) and (68) of [7] that cancelled each other. With 
the aforementioned conventions the correct identity is the one given by (93), (94), and 
Jy = + • ■ ■ (see Eq. (57) of [7] for the complete expression of Jy). 
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In order to go from (91) to (92) we look at the Ward identity (91) for the 
particular vector held 


V^ = K, 


{w,w) * 


(95) 


By the identity (93) the RHS becomes the right term in (92). That the LHS 


also coincides follows from an integration by parts together with 


a; -.). 


[z — w) 


(96) 


Similarly one can obtain the negative helicity soft theorem by looking at the 
Ward identity associated to the vector held 


To go from (92) to (91) one performs the operation —(dvr) ^ J d?wV^d% 


on both sides of (92). The LHS becomes the left term in (91) (for a vector 
held with vanishing component). The factors on the RHS are given by: 


(dvr)-! / (fwV^dlr^^ 

j {w,w) 


= -(d7r)-^T 


w 


with 




(97) 


(98) 


(here we used the fact that Jy depends linearly on the argument V, so that 
the operation f (PwV^d^ can be pushed inside the argument). Finally, using 


d-- 


wY 


w) 


= — z) 


(99) 


one hnds that W = V^dw. Thus one recovers the Ward identity for a vector 
held of the form V^du)- From the negative helicity CS theorem one can 
similarly obtain the Ward identity associated to a vector held of the form 
By adding the resulting positive and negative helicity identities one 


obtains the general form (91). 


6.3 External massive scalars 


Using the quantum version of (d6) 


[h.Q 


= iC 


Vh<- 


( 100 ) 
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the proposed Ward identity for external massive scalars takes the form: 


— lim (1 + EsOe, 
dvr Es^o+ 


ta, ta)5|in)+l/"'(?j^(out|a°'^*(ii^s, ta, ta)5|in)) = 


^i:;.^(out|5|in). (101) 


Here is the derivative V.^da on the i-th particle momentum variables 
(assumed to be outgoing; for incoming particles the factor is -Cy^). 

On the other hand the subleading soft theorem reads: 


lim (I+EsOe ) (out I a 


out 


(E„'«;,'«;)5|in) = ^ 


_ £+(w,w) - Pi + 


Q-Pi 


w, w)q^Jl"^ (out|5|in). 

( 102 ) 


with the total angular momentum on the i-th massive particle. We now 
repeat the steps that led to the equivalence in the massless case, starting 
in the direction soft theorem —)■ Ward identity. Accordingly, let us act with 
— (dvr)"^ J d?wV'^d% on both sides of (102) and with —(dvr)"^ J d'^wV'^d^ on 
the negative helicity soft theorem. Adding both terms we obtain an identity: 


— lim (1 + EsOe, 
dvr Es^o+ 


d w{V'^d:^{ont\a°^ {Eg, w, tc)5|in)+l/"'(?;;(out|a_ (h^^, tc, tc)5|in)) = 


^ £pi(out|5|in), (103) 


where 

H"(p/m) = J d‘^wG'^{p/m;w,w)V^{w,w), (lOd) 

with 

G%{p/m;w,w)da := ^ (105) 

and similar expression for the G° component. It can be verihed this function 
takes the following simple form m- 


G%{p/m;w,w)da 


3 VtM 

Stt {{q/Es) ■ {p/m)Y 


{q^lEg)dB{q''lEg)J^,. 


(106) 
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Finally, parametrizing the 3-momentum particle as p = mpx it can be shown 
G%{p,x;q) = G%{'p/m;q) satishes: 


{A-2)G%{p,x;q) = 0, 


Do,G%{p,x;q) = 0 , 


p^OO 


lim G^{p,x;q) = 6i6^^\x,q), 


(107) 


where A and act on the (p, x) coordinates!^ It then follows that G% is 
the Green’s function for Eq. (26), and so E = Vy,. Thus the identity (103) 
coincides with the Ward identity ( |101[ ). 

To go from (101) to (102) we repeat the steps as in the massless case: 
Consider the Ward identity for the particular vector held 


By the same argument as in the pure gravity case, the LHS becomes the left 
term in (102). The differential operator on the RHS is given by 

\2 


Vy = 


(Tw G7 {p/m]w,w) 


[w - z. 


{w - Zs) ’ 

e+{w,w)-p , _ 

= — - ^ , 


q-p 


(109) 

( 110 ) 


where we used the dehnition (105) of G% integrated by parts and used Eq. 
(96). One thus recovers the right term in (102). 


7 Remarks 


From the results of the present paper together with those of the Maxwell 
case m we arrive at the general picture that the ‘soft factors’ that multiply 
the hard particles in the soft theorems have the role of ‘potentials’ for the 
Green’s functions that implement the large gauge transformation on such 
hard particles (either massless and massive; in the massless the “Green’s 
function” is just an identity kernel). We are referring here to identities (64) 


and (72) for the case of supertanslations and identities (99) and (105) for the 


sphere diffeomorphisms (and similar identities for the Maxwell case [TT]). 
In each case, the 2-dimensional differential operator one is ‘inverting’ are 
the ones arising in the soft part of the charge associated to the large gauge 


lOn 


m —>■ 0 limit and using Eq. (99). 


The boundary condition in (107) can be obtained by casting the p —>■ oo limit as a 
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transformation nnder consideration. It is so far nnclear to ns what is the 
nnderlaying physical or geometrical mechanism behind snch identities. We 
hope to clarify it in fntnre investigations. 

Another generic feature of the Green’s functions and their potentials re¬ 
gards the relationship between the massive and massless cases: The iden¬ 
tity kernel “Green’s function” for the massless case can be obtained as a 
m —)■ 0, p = constant limit of the massive Green’s function. Mathematically 
this massless limit is equivalent to taking p —)■ cx) keeping m =constant. 
From this perspective, the fact that the Green’s function becomes an iden¬ 
tity kernel ensures the correct boundary value of the gauge parameter that 
is being extended from null to time-like inhnity. 

Finally, we point out a subtlety regarding these potentials: Whereas the 
Green’s functions are globally dehned on the q sphere, the potentials have a 
singularity at a point. The ‘potentials’ we used have the singularity at the 
south pole w,w = oo, which in turn is related to the fact that the expres¬ 
sion for the polarization tensor used is only valid for soft graviton directions 
away from this point. When showing a soft theorem from a Ward identity 
(both massless and massive case), an integration by parts is performed. The 
smearing function being in integrated are such that the would be ‘boundary’ 
contribution from the south pole vanishes. 

8 Conclusions and outlook 

For the case of pure gravity, due to the seminal work of Strominger and his 
collaborators and many works which have followed, a clear picture is emerging 
as regards to the symmetry of quantum gravity S matrix in perturbative 
regime and its remarkable connections with some well known factorization 
theorems in the quantum gravity literature. However all the work done so far 
was in the context of pure gravity or at most in the context of pure gravity 
coupled to massless matter helds. This was due to the fact that the classical 
symmetry group of asymptotically flat geometries is most clearly understood 
at null inhnity and is characterized in terms of so-called supertranslations 
and the Lorentz group (or its extensions like the local conformal symmetry 
or smooth diffeomorphisms of 3“^) acting on the conformal sphere. 

A candidate symmetry group of semi-classical quantum gravity S matrix 
when the scattering states are massless is Q and the associated Ward identi¬ 
ties are equivalent to certain soft theorems. A natural question then emerges: 
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Can we extend the above ideas to situations where perturbative gravity cou¬ 
ples to massive particles? In this paper we argued that this generalized BMS 
group Q can be made to act not just on null inhnity but also on time-like 
inhnity. A key ingredient for such extension is a description of time-like inhn¬ 
ity as a constant curvature spacelike hyperboloid that parametrizes all unit 
time-like vectors at inhnity. 

We showed that there is a way to map generalized BMS vector helds and 
obtain a group of large diheomorphisms which acts at the aforementioned no¬ 
tion of time-like inhnity and proposed this group ^ as a candidate symmetry 
group of the S-matrix. We showed that this group has a well dehned action 
on the free states of the theory and in fact the corresponding Ward identities 
are equivalent to Weinberg as well as Cachazo-Strominger soft theorems. As 
emphasized in the previous section, these equivalences rely upon intricate re¬ 
lationships between leading and sub-leading soft factors and Green functions 
which map supertranslation generators and sphere vector helds to functions 
and vector helds at the asymptotic hyperboloid respectively. 

A number of important issues remain to be understood. Here as in all 
previous papers we have worked with a formal S-matrix dehned by assuming 
that asymptotic states are free states of the theory. However as has been well 
known for a long time, in the presence of long-range interactions asymptotic 
states of the theory are not free but are obtained by dressing free state with 
a coherent cloud of soft particles. This tension is already visible at the 
classical level. If we analyze the equation of motion of a massive scalar held 
coupled to linearized gravity, its solutions in the asymptotic (r —)■ cx)) diher 
from the asymptotic limits of the free held by a phase factor exp(iyhrTlnr) 
where hir is the coefficient of ^ in a expansion of the linearized metric. 
The action of Q vector helds will in general be ill-dehned on such states. 
This is most clearly seen for the sphere vector helds which take the form 
Their action on dressed helds contains a multiplicative 
term proportional to \n{T) dah^rr ■ This implies that the action on the 
symplectic modes of massive helds, if we consider the helds to be dressed in 
contrast to free, is divergent. We believe that such divergences are intricately 
tied to the IR divergences that have been analyzed in [18] and these relations 
and their treatment remains an important question that remains unexplored. 
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A Equations (25) and (26) 


In this appendix we derive equations (25) and (26), which are the equations 


satisfied by the generators of Q at time-like infinity. We first derive the 


differential equations in (25), (26) and then derive the boundary conditions. 


A.l Differential equations 

We start with the assumed expansion of the vector field at time-like infinity: 


c ( t , p , x ) = r(p,i) + o(T') 


(111) 

( 112 ) 


and compute (18). In hyperbolic coordinates, the wave operator takes the 
form: 

□ = -V^ + (113) 

where is the inverse of the metric in ?{. The nonzero Christoffel symbols 
are: 


r;. = rv, n, = T-‘y, r 


( 114 ) 


ap ~ ' "‘afii ^ Pt ~ ' '^j3i ^ 0'yi 

with T g^ th e Christoffel symbols of the covariant derivative Da on 'H. Ap¬ 
plying (1131) to (lull), (|112l) one gets 


□r = r-i2Zl„r+ r-2(Ar + 2Zl„e^^)“ + 3r)+0 (t-^) (115) 

= r-2(AC-2C)+0(r-2). (116) 

he vanishing of the leading terms corresponds to the first two equations in 


(26). To obtain the equation for we need further information. For this we 


now look at the divergence of the vector field: 

V' := + T-‘(3r + + 0 (t-=). 


( 117 ) 


Now, from equations (18) and (19), it follows that ijj satisfies = 0 with 
fall-off "0 = 0{r~^) at null infinity. This is precisely the fall-off behavior of 
regular massless scalar fields. Such fields are known to decay at time-like 


infinity as 0 = 0{t ^), see for instance [Ejl^ The 0(r°) term in (117) is 


'^We thank M. Reisenberger for help on this point. 
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already known to vanish from Eq. (115). The vanishing of the 0{t term 
gives us the missing information: 




(118) 


which combined with the vanishing of the 0{t term in (115) gives the 


hrst equation in (25) 


A.2 Boundary conditions 


We now write the vector held (22), (23) in (u, r) coordinates to impose the 


boundary condition (20) at null inhnity. The change of coordinates is given 
by: 

u = r('\/l + - p), r = pr. (119) 

We want to look at the r —>■ oo, m = const limit. In this limit, we have 

p —)■ cxD, r —)■ cxD, 'r/(2p) = const., (120) 

with 


r/(2p) =u + 0{r ^). 


Condition (20) on the sphere components tell us: 

= v^ + 0{r-^). 


( 121 ) 


( 122 ) 


Writing the condition = 0 together with (122) hxes the asymptotic 

form of to: 

= -pa + 0{l), (123) 

where 2q: is the 2-d divergence of The r component of the vector held 
has the asymptotic form: 

r = (B,r)e + {drr)C = Tp + + pC) + 0(r-‘), (124) 

and has to satisfy the condition 

f = -ar + 0(l). (125) 


From (123) the leading —ar term is recovered and we are left with the con¬ 
dition: 

pe" + ^(i)P = o(i). ( 126 ) 
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We finally discuss the condition on the u component. Using 


_ = yiT? - p = (2p)-‘+ 0(p-=>), (127) 

I - + (128) 

one hnds: 

e = -r/(2p2)^P + ^7(2p) - e(')7(2p") + • • •, (129) 

where the dots indicate subleading terms. This needs to asymptote to: 


lim 7 = / + 

r^oo 


(130) 


From (123) we see that the hrst term in (129) reproduces the ua term in 
(130). We are then left with the condition: 


lim ^V( 2 p) - ^0)7(2p=) = /. 

P^OO 


(131) 


In order for (126) and (131) to be simultaneously satished, we must have: 


r = p/ + o(i) 

jdV = _p2/ + 0(p). 


(132) 

(133) 


To summarize, the asymptotic condition (20) translate into conditions ( 122[ ), 
(123), (132) and (133) for the vector held in (p, r) coordinates. 


B Equations (64) and (66) 


The tensorial structure in {w,w) implies acts as in Eq. (43). Performing 
this operation one hnds: 

du,{'y'^'^du,{'yy,ins{z,z;w,w))) = 1 + {^iTi6^'^\w-z)+2n{w-z)du,S^‘^\w-z)), 

-L ~f“ Z 2/ 

(134) 

where we used d^iw — z) ^ = 2'k5‘^{w — z). Seen as a distribution we have 
that {w — z)dw8^‘^\w — z) = —5^'^\w — z). Furthermore by the Dirac deltas we 
can set w = z in the prefactors of (134). It then follows that as a distribution: 




( 135 ) 






















which corresponds to Eq. (64). 

We now show Eq. (66). In the {z,z) variables s{z, z]w,w) is a scalar. 
Hence the second derivative acts as: — Tl^dz = Y^dzlzzdz- Per¬ 

forming this operation one finds: 

'^^"'dz{r{zzdz^{z,z\w,w)) = ^-^^;^^{4n6^^\w-z)-2n{w-z)dz6^‘^\w-z)). 

1 -I- 2;;^ ^ 

(136) 

Identifying {w — z)dzS^'^\w — z) = — z) one finds 

7 ''''dz(jzzdzs(z,z;w,w)) = 27r6^^\w - z). (137) 

which corresponds to Eq. 
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